Prictional drag in dilute bilayer 2D hole systems 
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We develop a theory for frictional drag between two 2D hole layers in a dilute bilayer GaAs hole 
system, including effects of hole-hole and hole-phonon interactions. Our calculations suggest signif- 
icant enhancement of hole drag transresistivity over the corresponding electron drag results. This 
enhancement originates from the exchange induced renormalization of the single layer compress- 
ibility and the strong dependence of single layer conductivity on density. We also address the 
effect of hole-phonon interaction on the drag temperature dependence. Our calculated results are 
in reasonable quantitative agreement with recent experimental observations. 
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Frictional drag measurements of transresistivity in 
modulation doped GaAs electron bilayer systems have 
led to significant advances in our understanding of den- 
sity and temperature dependence of electron-electron and 
electron-phonon interactions in 2D systems 0. Recent 
interest [^-^] has focussed on the role of electron corre- 
lation effects on the drag resistivity, which should vary 
in a systematic manner as a function of electron den- 
sity and temperature. In particular a recent experiment 
by Pillarisetty et al. [^) reports drag measurements in 
very low density and high quality hole bilayers, where 
Coulomb interaction (i.e., correlation) effects should be 
strong by virtue of the large GaAs hole effective mass 
(m,^ « 0.4m for holes compared with m* w 0.07m for 
electrons where m is the vacuum electron mass) and 
very low hole density (hole density p = 10 10 — 10 11 cm -2 
in Ref. ||], whereas typical electron densities used bi- 
layer drag measurements have been comparatively high, 
> 10 n cm~ 2 ). In terms of the dimensionless interaction 
strength parameter r s = (p7r) _1 / 2 m*e 2 /(?i 2 K) where k is 
the background dielectric constant, which measures the 
ratio of the potential energy to the kinetic energy in the 
interacting hole system, Ref. || explores the strongly cor- 
related regime of r s « 20 — 40 whereas the earlier electron 
drag experiments explored the weak coupling regime of 
r s < 3. 

The experimental findings of Ref. Q were striking, and 
motivated our work. They are: a. the low density hole 
drag is 2 — 3 orders of magnitude larger than the corre- 
sponding electron drag results published in the literature 
fj]~i|)0; b. There are some small (but systematic) devi- 
ations of the observed low density hole drag resistivity 
from the expected p ~ T 2 Fermi liquid behavior. The low 
density data of Ref. || seem to better fit a T 25 behav- 
ior at low temperatures; c. The observed p/T 2 behavior 
m Ref. §, plotted as a function of T for various bilayer 
hole densities, is qualitatively similar to the correspond- 
ing electron drag results in the sense that p/T 2 at a fixed 
density shows a peak at some temperature T p which de- 



creases with decreasing density - the peak in p/T 2 as a 
function of T at the lowest hole density p = 10 10 cm~ 2 is 
very sharp; d. For bilayers with unequal hole densities, 
the drag resistivity at a fixed temperature plotted as a 
function of the density ratio p\/p 2 decreases monotoni- 
cally and does not exhibit a peak at the balance point as 
it does in the corresponding electron case. This peak is 
believed to arise from the 2kp phonon scattering. 

These experimental findings become particularly in- 
teresting due to their possible relation to the collection 
of transport anomalies in 2D systems referred to as the 
2D metal-insulator transition (2D MIT) phenomena (The 
samples of Ref. Q exhibit 2D MIT in each layer at 
p k 8.5 x 10 9 cm~ 2 ). These anomalies, observed largely in 
the large r s regime, have raised doubts regarding the ap- 
plicability of Fermi liquid theory to two dimensional sys- 
tems of charges (electrons or holes) at the large r s regime, 
and this question has been widely debated in the litera- 
ture (see Q for a review). In that context, then, it is par- 
ticularly important to examine whether drag measure- 
ments in the large r s regime may be understood within 
the Fermi liquid framework. 

In this paper we attempt, and largely succeed, to inter- 
pret the hole drag data of Ref. Q within a Fermi liquid 
approach. The important inputs to our theory are a per- 
turbative expression to the drag resistivity, based on a 
perturbative treatement of inter-layer hole-hole scatter- 
ing rate, a Hubbard approximation for the single layer 
polarization operator, and the experimentally measured 
density dependence of the single layer conductivity. In- 
puts of more minor significance are form factors that ac- 
count for the finite thickness of the two layers and hole- 
phonon interaction parameters. 

Our starting point for the calculation of the drag re- 
sistivity is the following expression: 

(3 da i da 2 f q 2 d 2 q dw F 1 (q,Lu)F 2 (q,Ld) 
P ^a 1 a 2 djndmJ (2tt) 2 2^ sinh 2 (/3w/2) ' 

where Fi(q,uj) = lmIlu(q,LL))\ul2(q,uj)\. In Eq. (1), 



1 



[3 = 1/T is the inverse temperature (we use units such 
that ks = h = 2e = 1, except in final formulas); pi t 2 
are the hole densities in layers 1 and 2; oi^ are the con- 
ductivities of each layer; q is the 2D wave vector in the 
layer; nii/I±22 are the irreducible hole polarizabilities in 
each layer; and uf| is the dynamically screened effective 
interlayer interaction. Several comments are in place re- 
garding this expression: first, when the drag resistivity 
is derived from the Boltzmann equation ( [9|]l0|] , and see 
also the approximation ^ ss ^± [ s being made. 

This approximation is valid for well conducting layers, 
but becomes invalid at the low densities relevant here 
(see below). The need to replace j± by ^ is discussed 
in Ref. |0|. Second, the dynamically screened interlayer 
interaction ufE satisfies a matrix Dyson equation, 



u(g,w) = [1 - v t (q, uj)U(q, u)] v t (q,w). 
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We include in the bare interaction v t (q,uj) = v c (g) + 
v ph (q,u>) both the hole- hole direct Coulomb interaction 
and the phonon mediated interaction v p , which in- 
cludes the acoustic phonon propagator and the appro- 
priate hole-phonon interaction matrix element [p"3f . The 
phonon mediated interaction includes both deformation 
potential and piezoelectric couplings between the holes 
and the acoustic phonons using the standard hole-phonon 
interaction parameters for GaAs . Following the com- 
mon approximation we assume the inter-layer polariza- 
tion operators IIi2,Il2i to be zero. 

When Eq. (Q) is used to analyze experiments carried 
out on identical layers in the commonly explored regime, 
of high density and low temperature, it yields, 
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with £ being the Riemann zeta function. The assump- 
tions involved in getting from Eq. (Q) to Eq. (|j) are large 
inter-layer separation (Jc F d 3> f , qrFd ^> 1, with k F be- 
ing the Fermi wave vector and qxp being the Thomas 
Fermi screening wave vector), Drude relation between 
conductivity and density (cr i oc p,-), and a Random Phase 
Approximation (RPA) in which fl^ is replaced by its 
value for non-interacting electrons. The contribution of 
electron-phonon interaction is negligible in this regime. 
Eq. (||) is typically smaller than the experimentally mea- 
sured value, by a factor of 2 — 5. For the experiment of 
Ref. ||, the disagreement is much larger, getting as large 
as a factor of 500. 

Our calculation differs from that leading to Eq. (||) in 
five points, all of them leading to an increase of the drag 
resistivity. 

First, we use the Hubbard approximation to calculate 
the polarization operators fl^. Within this approxima- 
tion, 

Un{q,u) = [l + <(g)n°(q, W )G(g)]- 1 n^( g , W ), (4) 



where n°(q,w) is the polarization operator for non- 
interacting electrons, and G(q) = q/2^/q 2 + k F is the lo- 
cal field correction. For small q the Hubbard approxima- 
tion amounts to the introduction of a Landau parameter 
fo = —vf(q)G(q) w making the inverse compressibil- 
ity of each layer |^ = . This approximation is 

then consistent with measurements of the compressibility 
that yield negative values. In general, the effect of the 
Hubbard approximation is to decrease the bare intralayer 
interaction by a factor of I — G(q). The fact that the 
exchange-driven local field correction, G(q), affects only 
the intra-layer interaction is what makes its effect signif- 
icant: the interaction potential between anti-symmetric 
charge densities, vf 1 (q)(l — G(q)) — vf 2 {q), is modified 
from 2ne 2 d/n to ^-(d — -^-)- This increases the drag 
resistivity, and in the limit of k F d 1, q F Fd 3> 1, Eq. 
(p|) is replaced by 



PD 



£(3)7r h f k^Tn 



32 



? 2 k%d 2 



2k F d 
2k F d- 1 



(5) 



We note by passing that even for densities that are not 
very low, the fourth power of the last term in Eq. (|^) 
makes the correction quite significant. In the low den- 
sity regime we consider, the approximation Eq. (|5|) is 
not valid, since scattreing events with large momentum 
transfer have a significant contribution. Our calculations 
below show that in that regime the local field correction 
increases the drag resistivity by about an order of mag- 
nitude. 

Second, we do not use the approximation oi cx pi. 
Rather, we use measured values of cr(p), provided to us by 
Pillarisetty et al. jlj] , to extract 4 21 for use in our calcu- 
lation of the drag resistivity. At trie lowest density mea- 
sured ^^-f 1 ^ — 10 — 20, so this correction increases the 
drag resistivity by about an order of magnitude. These 
measurements of cr{p) are limited, at this stage, to the 
range of densities between 0.72 — 2.15 • 10 10 cm -2 , and to 
the temperature range T < 0.5K . While this correction 
deacreases with increasing density and temperature, it is 
still between 3 — 5 at the edge of the measured range. 
Our calculation for the range where experimental values 

of (y^r^j are available is presented in Fig. 1. The 

other figures do not include this factor, for lack of avail- 
ability of experimental data. 

Third, the assumption k F d » 1, implying that the 
relevant momentum exchange is much smaller than k Fl 
is not valid at the low density regime. Consequently, 
there is an important contribution due to large momen- 
tum transfer, hq ~ 2hk F scattering. Taking this contri- 
bution into account (by a numerical integration of Eq. 
(|l|)) leads to an increase of the drag by another factor of 
~ 2 as compared to the Boltzmann result. 

Fourth, the introduction of finite thickness form factors 
PH (see to the Coulomb and phonon interaction) 
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FIG. 1. The drag resistivity (a) and p{T)/T 2 (b) as a func- 
tion of temperature for various hole densities (p = 1.0, 1.5, 
2.0, 2.5 xl0 10 cm~ 2 , from top to bottom) calculated with all 
five correction factors as explained in the text. Throughout 
this paper we use a hole bilayer system with the layer sepa- 
ration of d = 300 A and the well width of a = 150 A 

effectively decreases d, and increases the drag resistivity 
by about a factor 2. 

And fifth, we include the phonon contribution to the 
interaction in addition to the Coulomb one. This changes 
the result by less then 50%. 

Combining all these five factors, we are able to ac- 
count for most of the results of the measurements. We 
account for the very large increase of drag, as compared 
to measurements of electronic systems. Our Fig. 1 is in 
good quantitative agreement with the measured data, to 
within a factor of 2. This type of agreement is similar to 
what is obtained in the small r s limit. Our analysis yields 
a leading quadratic temperature dependence of the drag 
in the limit T — > (at least as long as the conductivity 
a is temperature independent in that limit). However, 
our numerical integration of Eq. |l|, as presented in Fig. 
2, indicates that even at the lowest measured temper- 
atures the drag resistivity does not follow a T 2 depen- 
dence. There are a number of reasons for that. First, at 
the low densities used in Ref. || [shown as the inset in 
our Fig. 2 to be compared with inset (b) in Fig. 3 of Ref. 
||] the Fermi energy (1 — 2K ) is not much larger than 
the measurement temperature range. Second, even well 
below IK the phonon contribution to drag in the low 
density hole bilayers is quite substantial (in contrast to 
electron systems where the phonon contribution is typi- 
cally a factor of 10 3 smaller for small layer separations of 
d = 300A or so used in Ref. || ) . As such we believe that 
the experimental departure from the T 2 behavior 
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FIG. 2. The calculated total hole frictional drag p(T)/T 2 
as a function of T for various hole densities (p — 1.0, 1.5, 
2.5, 3.5, 5.0, 7.0 x 10 10 cm~ 2 , from top to bottom) within (a) 
RPA dynamical screening theory, and (b) HA. In insets we 
show the variation in the calculated peak temperature T p as 
a function of the hole density p. 

reported in Ref. || is essentially a manifestation of the 
fact that phonon effects remain significant in the exper- 
iments, and the asymptotic T 2 regime is hard to reach 
in hole systems. We find that our calculated p(T) at low 
temperature is well approximated by a T 2 4 behavior for 
p = 2.0 x 10 10 cm~ 2 and the exponent increases as the 
hole density decreases. In Fig. 4 we show our calculated 
contributions to the hole drag resistivity from individ- 
ual hole-hole and hole-phonon interactions as compared 
with the corresponding electron case. The importance of 
phonon effects to the hole drag transresistivity is mani- 
festly evident in Fig. 4. 

In Fig. 3 we qualitatively "explain" the particularly 
anomalous feature of the data in Ref. |(J [shown as the 
inset in our Fig. 2 to be compared with inset (b) in Fig. 
3 of Ref. ||], i.e., the non-existence of a peak in p/T 2 
as a function of the density ratio p\/p2- In qualitative 
agreement with Ref. || the calculated drag resistivity at 
a fixed temperature decreases monotonically as a func- 
tion of the hole density ratio p\ /p2 without showing any 
phonon-induced peak at the balance point p\ = P2, as 
has been observed in corresponding electron bilayer ex- 
periments §. This peak arises from the sharp 2D Fermi 
surfaces in the two electron layers which, when perfectly 
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FIG. 3. The drag resistivity in RPA for various density 
ratio (P1/P2, where pi = Pdrive and p2 = pdrag) with the fixed 
drag layer density (p2 = 2.0xl0 10 cm -2 ). Inset shows the drag 
resistivity as a function of a density ratio at T = 1.0K. 

matched at the balance point p\ = P2, leads to en- 
hanced phonon scattering, leading to the peak resistivity 
at pi = p2 ■ No such peak exists in the hole bilayer case 
because of the small Fermi temperature in the hole case, 
Tp = IAK, and the large Bloch-Gruneisen temperature, 
T BG = 2.8K, (for p = 2x 10 10 cm~ 2 ). Therefore typical 
T/Tp is rather large in the hole case, leading to thermally 
broadened Fermi surfaces in the low density hole bilayers 
of Ref. [|| which cannot exhibit any sharp Fermi sur- 
face effects, and consequently the so-called ll 2kp phonon 
peak" arising from the Fermi surface matching is ab- 
sent. The peak in pd/T 2 , appearing in Fig. 3, takes 
place at high temperatures, above the Fermi energy, and 
is unlikely to be related to the experimental peak that ap- 
pears around 0.5K. We believe it may be related to the 
two opposing temperature dependences of the hole-hole 
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FIG. 4. (a) The RPA-calculated electron drag re- 
sistivities from electron-electron v c (dot-dashed), elec- 
tron-phonon v ph (dashed) , and total interaction vt (solid) for 
n — 2. Ox 10 10 cm -2 . (b) The same for the hole drag resistivity. 



scattering rate on one hand and the (^^f-^j factor on 
the other hand. 

In conclusion we have developed a theory for frictional 
interlayer drag resistivity in low density hole bilayers in- 
cluding both hole-hole and hole-phonon interaction ef- 
fects. We explain the observed [|| dramatic deviation 
of drag resistivity from an RPA-based Boltzmann equa- 
tion calculation as arising from several factors, out of 
which the most important ones are the dependence of 
the single-layer resistivity on density and the exchange 
induced renormalization of the single-layer polarization 
operator. Our calculation explains the large increase 
in the magnitude of drag and the absence of phonon- 
induced 2kp peak in the drag between matched densities. 
We also obtain reasonable qualitative understaning of the 
temperature dependence of the drag resistivity. Our cal- 
culation are all done within a perturbative Fermi-liquid 
based approach. 
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